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We discuss a metal-insulator transition caused by random 
couplings of magnetic moments in itinerant systems. An ana- 
lytic solution for the single particle Green function is derived 
from dynamical self consistency equations, the corresponding 
density of states is characterized by the opening of a gap. The 
scaling behavior of observables is analyzed in the framework 
of a scaling theory and different crossover lines are identi- 
fied. A fluctuation expansion around the mean field solution 
accounts for both interaction and localization effects in a con- 
sistent manner and is argued to be relevant for the description 
of the recently discovered metal-insulator transition in 
2d electronic systems. 

PACS numbers: 71.20.-b, 71.55.Jv, 75.10.Nr 



Interaction driven metal-insulator transitions (MITs) 
as well as magnetism are among the most striking con- 
sequences of strong electronic correlations. For many 
years progress in the theoretical investigation of these 
phenomena was impeded by their manifestly nonpertur- 
bative nature, and only the introduction of powerful new 
methods like self-consistent dynamical mean field theo- 
ries (DMFT) has initiated progress ||,§. The modifi- 
cation of the disorder induced localization transition by 
interactions has been the subject of intensive studies ||, 
whereas the profound effects of random interactions on 
the band structure have attracted interest only recently 
H[7| in the form of a hard gap causing a crossover from 
variable range hopping to activated behavior in the con- 
ductivity. In this letter we explore the phase diagram 
of metallic spin glasses in the vicinity of the quantum 
phase transition between a gapped insulating phase and a 
metallic phase (SG-MIT) with a density of states (DoS) 
at the Fermi level increasing continuously from zero as 
the magnetic interaction strength is reduced. Applica- 
tions to the recently discovered MIT in 2d systems [02] 
and the spin glass phase in High T c superconductors E] 
are discussed. 

While in clean systems the MIT so far can only be de- 
scribed with the help of numerical methods, the preva- 
lence of disorder correlations in spin glasses makes it pos- 
sible to find an analytic solution of dynamical mean field 
equations. A Hamiltonian with hopping elements out of 
the Gaussian orthogonal ensemble leads in infinite space 
dimensions to a set of equations equivalent to the lo- 
cal impurity self-consistent approximation (LISA) Q as 
was already pointed out in . For the typical magnetic 
interaction J much smaller than the kinetic energy the 



DoS is semicircular with a width 2Eq and the system is a 
metallic spin liquid, whereas for Eq = (32/3n)J a quan- 
tum transition from paramagnet to spin glass takes place 
[p|-|ri| . A further decrease of Eq/J leads to a suppression 
of the DoS at the Fermi level pp until at Eq c = 2.59 J a 
gap opens up. The dynamical selfconsistent theory treats 
both the disordered electron hopping and the spin-spin 
interaction in a consistent manner and is therefor a good 
starting point for a fluctuation expansion. From the ap- 
pearance of replica symmetry breaking in the fluctuation 
theory it is conjectured that this transition is preceded by 
a Griffith phase in analogy with the ferromagnetic phase 
in the random field Ising model |12]| . Thermal fluctua- 
tions affect the SG-MIT in a double way: they smear out 
the transition itself and they shift its location by reduc- 
ing the magnetic correlations. These effects are described 
by crossover lines separating quantum critical, insulating 
and metallic behavior of DoS and transport coefficients. 
While in infinite space dimensions the SG-MIT has mag- 
netic correlations at its origin, in finite space dimensions 
localization fluctuations become important and the tran- 
sition has to be compared with the Anderson Mott MIT 
II . Both the existence of the order parameter DOS and 
the quartic critical theory Jl3| with a dangerously irrele- 
vant variable suggest that the similarity persists in fluc- 
tuation theory. However, we expect that the detailed 
treatment of interaction effects in our mean field theory 
and the correct choice of the ground state will render the 
renormalization group (RG) theory of the SG-MIT free 
of instabilities present in the usual sigma model approach 
||l4| , [l5f to interacting disordered systems. 
We study a single band model with the Hamiltonian 
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comprised of a kinetic part with real symmetric hop- 
ping matrix elements distributed according to a Gaus- 
sian probability measure with moments < ty >= and 
< (tij) >= M(i — j) , and a random interaction be- 
tween Ising spins with variance J 2 (SK model). It can 
be viewed as effective t—J model describing the inter- 
action of spin fluctuations at randomly distributed im- 
purity sites with strong Coulomb repulsion Jig ]. The 
disorder average is performed by means of the replica 
method, the replicated partition function is represented 
with the help of Grassmann integrals. The four fermion 
and four spin terms are decoupled with quaternionic ma- 
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trix fields R ab (x;T, r') and scalar fields Q ab (x; r, r'). In 
infinite space dimensions dynamical saddle point equa- 
tions become exact and the stationary values can be iden- 
tified with the single particle Green function and the spin 
autocorrelation function, respectively. After further de- 
coupling of the spin-spin interaction the fermions can be 
integrated out. For a replica symmetric and static ap- 
proximation of the linear susceptibility and in the limit 
of zero temperature two of the three integrals over spin 
decoupling fields can be solved exactly by the method 
of steepest descent. One ends up with a spectral type 
self-consistency equation for the the thermal single par- 
ticle Green function G(ei) — (A/iEfyn with n denoting 
the saddle point value of the charge decoupling fields. 
The shortcomings of this calculation are remedied in the 
framework of a Landau Ginzburg theory for the SG- 
MIT, which allows to include the effects of finite temper- 
ature, inelastic interactions and replica symmetry break- 
ing. The self-consistency equation for the Green function 
reads 
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where the weight function is given by 

A{z) = J-e[\z\- X0 \<t>'{z)\]{l - X *<t>"{z)]e-^- XaW{z) 
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The Edwards Anderson order parameter q = lim^oo < 
S(t)S(0) > and the zero frequency component \o of the 
local magnetic susceptibility have to be determined by 
extremizing the free energy, 4>(z) = Trlog[l + z 2 /(ri + 
e;) 2 ] is a functional of the electron Green function, all 
energies and magnetic fields are measured in units of the 
average spin coupling J. The derivation of eqs.(|]j|) wm 
be presented elsewhere. 

In order to gain insight into the behavior of the spectral 
type function A(z) it is useful to calculate it non self- 
consistently for a semi-elliptic Green function and dif- 
ferent values of the bandwidth Eq. Large values of Eq 
yield a functional 4>(z) with a maximum at z = and 
A(z) — Aq + O(z) for small z, consequently solutions of 
(|J) are metallic. For Eq < 1.5 however, <fi develops a 
double peak structure and therefor the weight function 
A(z) is gapped around z=0. In this case the denominator 
in eq.(^) can be expanded in powers of n + Keeping 
only the leading frequency dependence and adding the 
dominant temperature correction the result reads 
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and <5i is a positive constant. All expectation values are 
taken with respect to the weight function A(z). For the 
special value of the bare band width E = 2/ y/ < 1/z 2 > 
corresponding to S = one finds at zero temperature a 
solution of the form r; ~ sgn(ei) \/|ez| with a dynamical 
critical exponent z = 3. In order to locate the critical 
point we have performed a selfconsistent calculation of 
the Green function. To find an analytic expression for the 
functional $ the self-energy is modeled by ri = k(^) 1 ^ 3 
for \ei\ < A with a cutoff A determined by the normaliza- 
tion condition for the DoS 1 = J_ A p(e). A free propaga- 
tor G(ei) = 4- is used for energies larger than the cutoff 
A. The functional $ is in this approximation given by 
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Using this functional we have solved numerically the sys- 
tem of coupled selfconsistency equations for critical band 
width, local susceptibility, and spin glass order parame- 
ter, and found the result E ac — 2.53J, \oc — 0.79/ J, 
and q c = 0.66. We have verified that a solution ~ ^JTi 
corresponding to a parabolic gap in A(z) does not satisfy 
eq®. 

To discuss various crossover scenarios one needs the full 
solution for the single particle Green function, though. 
It is given by 
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FIG. 1. Single particle density of states DoS versus energy 
for different values of the electron hopping strength corre- 
sponding to metallic, critical and insulating behavior. 
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Here S — 5q — TSi plays role of a mass squared which 
separates the disordered (insulating) phase with pos- 
itive 6 and vanishing pp from the ordered (metallic) 
phase. The dependence of the DoS on the bare band- 
width Eq is displayed in fig.(l). The temperature de- 
pendence of 8 describes the effects originating in the re- 
duction of magnetic correlations by thermal fluctuations. 
The DoS at the Fermi level depends on temperature via 

Pf = — %z\ — • For low T the constant term Sn dom- 
inates and the DoS is approximately constant, whereas 
for T ^ thermal fluctuations take over and the sys- 
tem behaves as if it were at its quantum critical point. 
Hence for all observables dominated by the critical DoS, 
the quantum critical (QC) region is delimited by the lines 
T = JM S ee fig.(2). _ 

Given the fact that in the replica symmetric solution of 
the SK-model (as well as in simulations of finite dimen- 
sional spin glasses) the distribution of local fields P(hi oc ) 
is nonzero for hi oc = it is questionable whether the 
frozen in magnetic moments in a spin glass can cause 
a gap in the single particle DoS and drive the system 
insulating. The investigation of generalized Thouless- 
Anderson-Palmer (TAP) equations JtJJTtJ suggests that 
a spin at site i lowers its energy by polarizing its neigh- 
borhood. The resulting self energy E g = [J2j Xjj(Jij) 2 }av 
favors the single occupation of sites and thus insulating 
behavior. The influence of replica symmetry breaking 
and finite range interactions on the magnetic self energy 
are discussed in [7|, [l^ , [T^] . 

For finite dimensional systems the DC conductivity is 
calculated to leading order in a 1 /^-expansion |2Cj from 
the density-density response function D(q,uj n ) utilizing 
the Kubo relation 

d 

udc — lira, lj \\m. -—rImD R (q,uj) . (8) 

The result o\dc = e 2 pfD with the charge diffusion con- 
stant D = 2irb m pF (6m being the second Fourier coeffi- 
cient of the hopping correlation function M(q)) exhibits 
a linear dependence on T in the QC region since there 
the DoS is suppressed like VT. 

The results of the last paragraphs can be collected in the 
unifying framework of a scaling theory In the finite tem- 
perature formalism the DoS p(ei) — ~^§^ is the deriva- 
tive of the free energy with respect to an external fre- 
quency e;. This external frequency is analogous to an ex- 
ternal field in the context of a magnetic phase transition, 
a nonvanishing pf breaks the symmetry between ad- 
vanced and retarded energies. The transition is sharp for 
zero external field, whereas in the general case the DoS 
obeys the scaling relation p ~ \5\^1Z±{j^^) . The limit- 
ing behavior of the scaling function is \vcn x ^ 00 'R.±{x) = 
x x / z and lim x ^o TZ-(x) ~ #(t1)- The two relevant scaling 
fields T and Eq — Eq c are related by a crossover exponent 
<t> = 1 and hence can be combined in a single variable 



5. From the denominator uo n + Dq 2 of the diffusive two 
particle propagator one determines the scaling behavior 
of the diffusion constant as 

D~<T(*- 2 ) , (9) 

and inferring the violation of hyper-scaling derived from 
the RG analysis we finally obtain 

So far all results have been obtained in the framework 
of a replica symmetric and spin static dynamical mean 
field calculation. In the following we will argue that the 
physics of a finite dimensional system with replica sym- 
metry breaking and inelastic many-particle interactions 
is correctly described by a dynamical Landau-Ginzburg 
theory. 
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FIG. 2. Phase diagram in the vicinity of the SG-MIT at 
a bare bandwidth Eo = 2.53J. The transition lines separate 
the quantum critical region with a T 1,/2 -dependence of the 
DoS from the insulating regime with a gap at lower Eo and 
from the metallic region. 

In order to calculate critical properties of the finite di- 
mensional SG-MIT we have expanded the free energy in 
fluctuations about the dynamic mean field solution. The 
fluctuation fields are hermitian An x An matrices which 
depend on two time and one space variable. They obey 
charge conjugation symmetry in the form CR T C T = itf 
with the charge conjugation matrix C_ = ia-2 ® cto- We 
obtained for the quantum action governing the critical 
theory after a suitable rescaling of space and time vari- 
ables 

A=\f d d x{T 2 £ R&to e k , £J )[(-V 2 + e 2 '* + ^ 

J ei,e k 

+el /3 e¥ 3 + l)5 aS S^ - 7 ^7|, 7 ]< 5 (z; e h e k )} 

+uT 3 KUx; e k , ei)rf, a rf s R^{x; e m , e k - e t + e m ) 

ek,ei,e m 
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Here "f z = o~o®o~ z is a spin matrix and r ab depends on the 
Parisi solution of the spin glass problem. The An 2 critical 
modes satisfy Bpj z = y z R. Note that the quantum inter- 
action u has no influence on the mean-field solution as for 
a spin glass there is no net magnetization. The spin glass 
order parameter q ab modifies the random field term and 
describes the influence of frozen magnetic moments. The 
dangerously irrelevant variable t is already introduced at 
tree level to account for the shift of the upper critical 
dimension caused by disorder correlations. Due to the 
presence of replica symmetry breaking the value of its 
exponent 9 t will be different from 2 already in a two loop 
RG calculation. From the presence of replica symmetry 
breaking the appearance of a glassy phase preceding the 
SG-MIT can be conjectured in analogy to Q. 
In a tree level RG-analysis the quantum mechanical 
many-body interaction u is found to be irrelevant with 
scaling exponent — 9 U = —I in the vicinity of the upper 
critical dimension dc = 6. Actually u is dangerously ir- 
relevant and changes the crossover line between metallic 
and quantum critical region from the naively expected 
T ~ S zv to T ~ flC*")/^*"). As long as the quantum 
dynamics are irrelevant the lower critical dimension cfp 
of the model is expected to be equal to two. However, 
if the interaction term becomes relevant in low spatial 
dimensions as it happens in the Finkel'stein theory of 
interacting electrons, it will shift d$ below two. This 
effect of a reduced lower critical dimension is well known 
in models with relevant quantum dynamics. Recently a 
MIT has been observed in high- mobility metal-oxide- 
semiconductor field-effect transistors (MOSFETs) 
and a variety of other 2d systems. These experiments 
challenged the widely held belief that two is the lower 
critical dimension for MIT's in all types of disordered 
systems. The RG theory of the MIT in interacting sys- 
tems 14,n| indeed admits the possibility of a metallic 
phase but is hampered by the occurrence of a magnetic 
instability in the triplet channel. A transition to a fer- 
romagnetic state seems unlikely as the 2d metallic phase 
is quenched by a parallel magnetic field p2] |. Since ran- 
domness is relevant in these materials it is reasonable to 
expect spin glass freezing of magnetic moments. The di- 
vergence of the triplet amplitude in the RG theory [ ^4|Jl5| 
occurs at a finite length scale, therefor the frozen mo- 
ments are most likely not single spins but ferromagneti- 
cally ordered clusters. These ordered clusters may form 
in regions with strong disorder, interact via a RKKY type 
interaction mediated by conduction electrons, and form 
a Stoner glass . The freezing of magnetic moments as 
represented by the random field term in eq.(ll) reduces 
the magnetic susceptibility and is expected to control the 
runaway flow in a RG analysis. The validity of these ideas 
can be tested by performing magnetic measurements and 



looking for irreversibilities indicative of random freezing 
of magnetic moments. 

In perovskite materials with narrow bands and a pro- 
nounced spin glass phase like LaSrCuO on the other hand 
there should be a possibility to observe the interaction 
driven SG-MIT by lowering the temperature below the 
spin glass freezing temperature. 

In summary, we derived dynamical selfconsistency equa- 
tions describing itinerant systems with random mag- 
netic couplings. The analytic solution displays a metal- 
insulator transition due to magnetic correlations. The 
mean field theory served as a starting point for the deriva- 
tion of a fluctuation theory describing the combination 
of localization and interaction fluctuations. The effects 
of randomly frozen magnetic moments are described by 
a random field term, and from the appearance of replica 
symmetry breaking in this term the presence of a glassy 
phase preceding the MIT was conjectured. 
Acknowledgment. We thank C. DiCastro and V. Pudalov 
for valuable discussions. 



[2 
[3 

[4; 

[5. 
[6 

[<: 

[8 

[9 

[10 

[11 

[12 
[13 

[14 
[15 



[16. 

li- 
tis: 

[19 
[20 

[21 

[22 



S.V. Kravchenko, et al., Phys. Rev. B 51, 7038 (1995); 
Phys. Rev. Lett. 77, 4938 (1996). 

D. Simonian, et al., Phys. Rev. Lett. 79, 2304 (1997). 
D. Vollhardt in Correlated electron systems, World Sci- 
entific, 57 (1994). 

A. Georges, G. Kotliar, W. Krauth, and M. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

Belitz and Kirckpatrick, Rev. Mod. Phys. 66, 261 (1994). 
Peihua Dai, Youzhu Zhang, and MP. Sarachik, Phys. 
Rev. Lett. 69, 1804 (1992), and references therein. 
R. Oppermann, B. Rosenow, Phys. Rev. Lett. 80, 4767 
(1998). 

F.C. Chou, et al., Phys. Rev. Lett. 75, 2204 (1995). 

R. Oppermann and M. Binderberger, Ann. Physik 3, 494 

(1994) . 

S. Sachdev, N. Read, and R. Oppermann, Phys. Rev. 52, 
10286 (1995). 

A. Georges and A. Sengupta, Phys. Rev. B 52, 10295 

(1995) . 

M. Mezard and M. Monasson, 



:ond-mat/9406013 



T. R. Kirkpatrick and D. Belitz, Phys. Rev. Lett. 74, 
1178 (1995). 

A.M. Finkelstein, Zh. Eksp. Teor. Fiz, 84, 168 (1983); Z. 
Physik 56, 189 (1984). 



C. Castcllani, C. DiCastro, and P.A. Lee, cond 



mat/9801006 



J.A. Hertz, Phys. Rev. B 19, 47 96 (1979) 



M. Rehker and R. Oppermann, :ond-mat/9806092. 
R. Oppermann and B. Rosenow, Europhys. Lett. 41, 525 
(1998). 



R. Oppermann and B. Rosenow, xmd- mat/9803239 



F. Wegner and R. Oppermann, Z. Physik B 34, 327 
(1979). 

M. Mezard and A. P. Young, Europhys. Lett 18, 653 
(1992). 

V.M. Pudalov, et al., JETP Lett. 65, 932 (1997). 



4 



